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ABSTRACT 


Markov Chain models have been used to forecast stocks in 
a wide range cf manpower systems. Studies have been done in 
many areas such as education planning, hospital planning, 
manufacturing, private research and development, a women's 
military unit, the civilian work force supporting the U.S. 
Navy and a state police organization. This study looks at 
such systems under conditions of change and develops’ the 
equations that describe the steady state distribution of 
personnel. The conditions of change include systems where 
recruitment is constant, increasing (decreasing) additively, 
or increasing (decreasing) multiplicatively and systems 
where the changes in total system size are additively or 
multiplicatively increasing (decreasing). Numerical exan- 
ples utilizing these models are provided, along with a 
Computer program of the formulas written in the language 
APL. 
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I. BACKGROUND 


A. INTRODUCTION 


According to Bartholomew and Forbes [Ref. 1: p. 1], 


Manpower planning is often defined as the attempt to 
| edae the supply of people with the jobs available for 
em. 


This is a concept that is an integral fart of military force 
planning, today. All branches of the Department of Defense 
are involved with studies concerning manpower planning. 
Since large numbers of personnel are involved, the use of 
Manpower models to determine the results of changes to the 
existing system before they are actually implemented can 
prevent costly future problems. This is especially true in 
terms of having the rignot number of personnel in the rignt 
places with the right grade ievels and skills, at the right 


time. To qucte again rrom Bartholomew and Forbes [Ref. 1: 
P- 6], 


A manpower model is a mathematical description of how 
change takes place in the systen. 


Much work has been done in the area of manpower/ 
Organizational modeling. The text by Bartholomew and Forbes 
(ceferenced above and hereafter referred to as Bat) compiles 
information cn the subject from numerous disciplines and 
Standardizes notation to describe some of the basic models 
and statistical techniques used in manpower planning. A 
central theme of this work is the idea of depictiny an 
Organization as a system of stocks and flows [Ref. 1: pp. 


3,4]. ‘The personnel within the organization under study are 


divided into classes, or categories, by specific character- 
istics (for example: age, time in service, grade, etc). 
The numbers of personnel within categories at a specific 
time are considered stocks. Flows are the movements of 
personnel from one category to another during a unit 
interval of time, such as a year. 

Various mathematical models that predict future stocks 
from one year (or other time period) to the next are 
described in BaF. Of particular interest here 1s the tran- 
sition model that is based on the theory of Markov Chains. 
In addition to describing the model {Ref. 1: pp. 85-98], BaF 
provides a BASIC computer program [Ref. 1: pp. 248-260], to 
calculate future stock vectors. That program was converted 
into an APL rogram at the Naval Postgraduate School (NPS) 
and is available on the NPS mainframe computer, an IB” 
system 370. The APL program, like its predecessor, computes 
stocks from cne year to the next based on input such as, 
initial stocks, recruitment and wastage (or attrition) rates 
and the transition probabilities for the stocks. This can 
be done repetitively and theoretically could compute future 
stocks forever. 

There are two versions of the basic Markov Chain model 
that are well developed in BaF. One is based on a fixed 
flow of recruits into the system andthe otner on a fixed 
total system size. Also of interest are other versions (or 
Options) that deal with a system undergoing growth where 
either recruitment or total system size grows in an additive 
or multiplicative fashion. The APL program allows the 
computation of future stocks for the cases where: 

(1) cecruitment is fixed, grows additively or muiti- 
plicatively; 
(2) the total system size grows additively or muiti- 


plicatively. 


Steady state is the condition which a system reaches, 
after sufficient time has fassed, so that the stocks no 
lorger change. There have been steady state models devel- 
oped for the cases where there is a fixed flow of recruits 
into the system and where the total system size is fixed. 
In both these cases the number of perscnnel in each category 
remains fixed after steady state 1s reached by the systen. 
The vector whose components are these numbers 1s called the 
steady state stock vector (S.S.S.V.). 

Under conditions of growth, the steady state of a systen 
must be redefined. In those cases, steady state does not 
exist inthe above sense, because the system 1S undergoing 
constant expansion and so stocks will never reach a fixed 
Size. However, it turns out that the stock sizes do reach 
constant propcertions and so it is possible to talk about a 
Steady state distribution vector (S.S.D.V.). Components of 
this vector represent the proportions of personnel that 
remain in each category forever once steady state has been 
reached. In effect, in taose cases where the distribution 
of personnel among the categories no longer changes with 
time the systean has reached a steady State oa 
distribution." 

Although, at the present time, the APL program at the 
Naval Postgraduate School predicts future stocks for a 
System that grows under conditions mentioned above, itis 
able to supply a steady state stock vector (S.S.5S.V.) only 
for the case where there is a fixed flow of recruits into 
the system. The only way to get the steady state distribdu- 
tion vector under the various conditions of growth 1S to 
continue to run the program, i.e., to forecast stocks for as 
Many years into the future as necessary until steady state 
is reached. This is not a Satisfactory method for several 


EcasOn'S > 


(1) it takes too long to run the model when the systen 
reaches steady state only after Many years’ have 
passed, 

(2) it is a trial and error process of going Dackeaae 
forth, in time, to try and find out at what point 
in time steady state has been reached, 

(3) occaSionally, even if the program 1S run out as 
tac as the present formatting of the prcgran 
allows, there can still be instances where steady 
state has not yet been reached. 

This thesis uses the ‘basic prediction equation't, as 
described in BaF [{Ref. 1: pp. 86-88], to correct the arore- 
mentioned deficiencies by developing models to compute 
analytically the steady state distributions for systens 
under varying conditicns of growth. Also developed was an 
APL progran, to pe added aS a sSutroutine to the fain 
program, that computes the S.S.D.V. for these options. 

The next section in this chapter sets out the necessary 
notation and equations to enable the reader to follow the 
development of these steady state models. For more amplifi- 
cation on the caSic Markov model see BaF [Ref. 13 pp 3-8, 
86-132]. 


Be. NOTATION AND BASIC EQUATIONS 


Stocks, representing personnel at an organization, are 


divided into k categories. The stock vector 


a(t) = [ n, (t) eho (t),- 2 = ¢ Ny (t)) 
iS a row vector where each element ny (t), represents tke 
humber of personnel in category i at time t. Time 1S meas- 


ured discretely with a time interval of unit lenyth cenoted 
(t= eee The most frequently used feriod of time is the 
year although it 1S Sometimes convenient, or necessary, to 


uSe quarters, months, etc. 


10 


Flows describe the movements of personnel. There are 
two basic types of flows in a system. The first type of 
flow is internal to the system and can be thought of as 
transfers from one category to another within the system. 
Demotions and promotions would be examples of this type of 
flow. These are denoted by the k by k matrix, N(t-1), where 
each element, Sy om uUe represents the number of personnel 
moving trrom category 1 to category j during the time 
interval (t-1,t). 

The other type of flow involves the two-way transfer 
between the system (organization) and the outside world. 
First, there are wastage (or attrition ) flows of the number 
of personnel leaving the system from category i, during the 
time interval (t-1,t), denoted by De cay (t-1) « (2 See? 6 ceeten K 
gmat = 1,2,...). Then there are recruitment (or accession) 
flows of personnel into category j from the outside world, 
during the period (t-1,t), denoted by No; (t). (7) Saale zane, K 
mame t = 1,2,<<)- Following BaF [Ref.j 1: pp. 3,4 j, the 
notation for recruitment flows is slightly different from 
that for internal fiows and wastage flows. Tiss tO 
remind the reader that in most organizations attrition and 
internal flows (transfers) are accounted for rfrirst and then 
recruitment is determined at the end of the time period. 

The Markov Chain model uses flow rates (or protabili- 
ties) rather than actual numters. P, a k by k matrix of 
probabilities has elements 

Pi eee Te : oe ae 

There are certain basSic asSumptions that a Markov Chain 
model needs tc meet to be valid: 

(1) the transition probabilitiés, Pe. do not change 


over time; 
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(2) the population in each category is homogeneous so 
that Pij represents the probability of each indi- 
vidual in category i moving, independently of any 
other individual, to category j. 

While these assumptions may not always be met, there 
have been enough studies done using the model to show they 
are sufficiently realistic to justify the use of the model. 
Even in cases where these assumptions don't hold very weil, 
the model was fcund to be quite useful for predicting future 
stocks [Refs. 2,3]. 

Next, W = (Wy Wo gooee W,) is defined as the wastage 
(attrition) rate vector where each element, w:, is aefined 


L 
as 


Become 
. ett 1) 


Ncte that, 
K 


C 


dP, i W. = 1 g 1 = dg. Zepearerere Ke o 
yet 


Finally, £ = (©, ,f>4---,f,) is defined as the requis 
ment proporticn vector where each element 


J R (t) 


is the proportion of total recruitment that will be allo- 


cated to category j, with 
K 


Ce) > No) (t) 


j=l 
being the tctal recruitment during the time interval 


(t=1,t). Lias Anpires Cer dis 
j=l 
BaF {Ref. 12 pp. 87,88], snow that future stocks may be 


calculated uSing the equation 


12 


K 
° (ct) = 2, ni top, + R(t)r ; (eqn 1.1) 
for j= 1,2,---,k and iC She This eguation states, 
that the stocks in category j at time t are equal to the 
number of personnel that moved into category j from anywhere 
within the system during tite period (t-1,t), plus’ the 
Number of personnel recruited into category j during’ the 


Same period. In matrix form the equation is 
n(t) = n(t-1)P + R(t)x See eles, F (Zen ay 


BaF (Ref. 1,: p- 88] refers to this as "the basic prediction 
eguation." . 

Another eguation dealing with changes in total systen 
Size 1s also developed in BaF [Ref. 1: pp. 94-96]. To see 


this, note that . 


N(t) = is (t) 


c=l 
is the total system size at time t and 


M(t) = N(t) - N(t-1) (eqn 1.3) 


represents the change in the total system size during the 
time interval (t-1,t). This implies that total recruitment 
during (t-1,t) can be expressed as: 
K 
R(t) = Da (ete, +98.({t) 2 0 ‘ (egn 1.4) 
Lf 
The first term in this expression 1S recruitment that 1s 
done to replace those who left (attrition from the various 
categories) while the second term consists of recruitmert 
midtewrs due to the change in total system size. 
Pioeitt) = 0, then only those personnel who Jleave the 


system are replaced, which implies that the total systen 
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size remains fixed. If M(t) > 0, new personnel are supplied 
to cover an increase in total system size, i.e., to fill new 
jobs. eee t ace then some jots must be eliminated as 
soon as they are vacated. Of course, R({t) 2 0 always. 
Referring back to Equation 1.1, replace R(t) using 


Eguation 1.4 and let 


Ge. = p.:. + “WE 
‘J 4 c J 
Then, as BaF [{Ref. 1:3 pp. 94,95], shows 
K 
ae = : By eid ela + (US) 3 ; (eqn 1.5) 
t=/ 
for j= 1,2,..-,k and t = 1,2,-.. . To put this 1HEORaaaE, 
nctation let 
= {q--} r Pp 3 = UB 6 ps P 


= 


Or 


Then Equation 1.5 may be expressed as 


De) = Gea) Ce ee, tS ie eee (eqn 7.6) 


Equations 1.2 and1.6 form the basis from which five 
Steady state sukmodels (or options) will be derived. The 
next section introduces these options and the direction of 


development. 


C. OBJECTIVES 


The objective of this thesis 18S tc develop steady state 
equations modeling five different. conditions of growth. 
Each of these conditions will be addressed separately as a 


“drOowtn. “Opt L.Cnw The first three options deal with 
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different types of recruitment. They use Equation 1.2 where 
R{t) 1s the total recruitment during the time period (t-1,t) 
and in these three options, R(t) is the only term changing. 
The last two cptions deal with changes in total system size. 
They utilize Equation 1.6 where. M(t) represents the change 
in total system size during the period (t-1,t) and in these 
two options, M(t) is the only term that changes. A complete 
list of these five options is: 
Give Erixed Lecruitment: R(t) = FP, t21, j&R>0O fixed, 
1.e., recruitment is always constant; 
(2) Additive increase (decrease) of recruitment: 
R(t) = R + (t-1)M, t21 
where Mis the constant amount of yearly increase 
(decrease) ; 
(3) Multiplicative increase (decrease) of recruitment: 
P(t) = Ret-1, t21, 
where 100 (e-1) is the percent of yearly increase 
(decrease) ; 
(4) Additive increase (decrease) in system size: 
M(t) = M , t21 
where Mis the constant amount of yearly increase 
(decrease) in total system size; 
(5) Multiplicative increase (decrease) in system size: 
M(t) = (6-1) 6t-1N(0) Pe ora, 
where 100 (6-1) is the percent of yearly increase 
(decrease) in total system size. 


i= Jeane option 2 or 6 = 1 1m Option 3, then those 


Options revert to option 1, of fixed recruitment. iijcaees — 0) 
mmoptten 4 Gf © = 1 ah option 5, then M(t) = 0 and the 
result iS a system with a fixed total system size. That 


case will be dealt with as subcases of options 4 and 5. The 
fixed total system size case and option one, where Ris 
fixed, are the only cases where it is possible to achieve a 


Steady state stock vector, i.e., where the stocks don'ft 
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change with time. In all other cases steady state will be 
considered to have been reached when a steady state distri- 
bution of personnel is achieved. 

Chapter Ii addresses, in depth, the analytical develop- 
ment of the steady state equation for each of the options. 
It starts with an account of the work done on option 1 in 
BaF. Then it proceeds with the derivation of the steady 
state equations for options 2 through 5 from either Equation 
1.2 O©” Equationeicce. It also presents exaaples, one a 
four-~grade hierarchical systen, the other a three grade 
non-hierarchical system. In addition to formulating analyt- 
ical models fcr options 2 through 5, an APL program has also 
been developed that will compute the steady state distribu- 
tion tor tnose options. This program is contained in 
Appendix F. Appendices A through E are computer printouts 
of terminal sessions of the two examples for each of options 


tT throughs. 
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Ii. OPTION DEVELOPMENT 

The options mentioned in Chapter I will be discussei 
here under the usual restrictions of manpower systems. 
These restrictions imply that w,> 0 for all i=1,2,...,k, 
since in manpewer systems attrition must be allowed from any 
category. This means that the P matrix, as defined in 
Chapter I, iS composed of non-negative elements and has [row 
sums all strictly less than one. The OQ matrix, also defined 
in Chapter I, 1S a Stochastic matrix (i.e., it is composed 
of non-negative elements and has row sums all egual to one). 
These facts along with the following two mathematical theo- 
rems are needed in the derivation of the steady state 


results for the five options. 


Theore 


——— = — 


4 


Ir P is a matrix composed of non-negative elements and 
row sums strictly less than one and@is a scalar then 
the matrix @I-P haS a unigue inverse for all vaiues of 
@ 21 fRef. 4: pe. 42 |}. 


heoren 2 


If Q is a stochastic matrix and @ 1s a scalar then the 
Matrix @I-Q has a unigue inverse for all values of 
8 > 1[Ref. 4: pp. 60,61 }. 


Ao RECRUITMENT OPTIONS 


The follcwing three options are all variations of 
Eguation 1.2, where k(t), representing total recruitment, 1S 
the term that takes on different forms. Conseguently, the 


Steady state equations for these options are also obtained 
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in somewhat different forms. These steady state eguations 


acre derived in the next three subsections. 
1. Option 1: Fixed Recruitment 


The steady state equation for this first option is 
developed in BaF {Ref. 1: p.90]. The analysis presented 
here is a brief summary of that in BaF, using the same 
assumptions. Recalling “the basic prediction equation, 


(Equation 1.2), and assuming that 


limit R(t) = R is fixed 
t— co 


it rollows that the steady state stock vector, 
n= limit = no(t) 
t—-co ~— 


also exists [Ref. 4: pp.40-43,48-50 J. Therefore, letting 


t— co in Equation 1.2 resuits in 

Or 

Pinally, knowing the matrix I-P has a unique inverse (see 
Theorem 1), the result | 


Be) S82 Bale 2) ha aie (eqn 251) 


is obtained for the steady state stock vector (S.S.S.V.).- 
After steady state has been reached, it follows that 


o(t) =n , 


for "large" t. The following two “examples will be used to 


illustrate the use of this steady state equation. 
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Example (a) 

This first example uses data from Example 4.6 in BaF 
{[Ref. 1: p. 97), which is taken from the women officer's 
system of one of the British services. ifwnets a hierbarch-— 


ical, four category systen. 
The inital stcck vector is: 
n(0) = (129, 74, 28, 11) 3 


The transition matrix 1S: 


.728 .102 0 0 
0 .83 .046 0 

-— 0 0 Bey 025 
0 


0 0 - 902 ; 
The total recruitment is: 
Hoe jo5, caixed for ali t>] . 
The recruitment proportion vector is: 


c= (ls 0, 0, 0) = 


Then 
27 2s LOZ 0 0 
0 me -~.046 0 
iIn-pP= 
0 0 Bl 33033 
ne 0 0 -098 
and the inverse is 
326 (ome. 2 0G 263 Spee | 
0 Secrewe. 26034 685 
Cit 
0 0 7.519 Zoe 
0 0 0 10.204/ . 


i 


Using Eguation 2.1 


3.6/6 2.206 - 763 5 a7 
0 Seco2 2-034 Ou 
0 0 Tsou 22552 
0 0 0 10.204 


(35) (1, 0, 0, 0) 


| 
il 


(128.66, 77521,926. (05,ec 3 7 ee 


After rounding to the nearest integer the 5.S.S.V. 1s: 
h = (129, 717527 : 


The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 
can be found in Appendix A, confirming the steady state 
stock vector above. As the printout shows, steady state is 


achieved by t = 10. 


Example (b) 


This 1S a non-hierarchical, three-grade systen. 
The inital stock vector is: 

n(0) = (300, 200, 1900)1 = 
The transition matrix is: 


S00 .venloiesss Os 
P = ae ee 0 
21S 0 «90 - 


The total recruitment is: 

Rit) =. 100, fixed for SEU td 
The recruitment proportion vector is: 

£ = (210, .048225)- 
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Then 


« 4 =<. 1> =—.05 
I-P = ~.75 a 0 
-.05 210 
and the inverse 1S 
4.267 18 22133 
(I ~- P)y-i = 4 2 2 
2uN433 4 11.066 ‘ 
Using Equation 2.1 
4.267 2.133 
ma(100)(.70, .07, +23) 4 2 2 


(mes peige2, 417-8) « 


After rounding to the nearest integer the S.S.S.V. is: 
mee (376, 79, 418) ™ 


The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 
can be found in Appendix A, confirming the steady state 
stock vector above. As the printout shows, steady state is 


achieved by t = 50. 
2- Option 2: Additiv 


This option addresses the fixed additive increase 
(decrease) in recruitment for each time period. Here, 


R + (t-1)M, ceplaces R(t) in Egquation 1.2 resulting 12 


n(t) =n(t-1)P + Ro + (t-1)Mc - (eqa 2.2) 
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= 


In this option, M 


would be negative recruitment (R(t) <0) 


sense in manpewer nodels. 


cases. 


Case 1: 


linit 
t —> 0 


which is identical to option 


Case 2: M > 0 


Assume that, 








re n(t) 
eet 
eXuS ess This is the same 
n (t) behaves like tMn . 
Equation 2.2 can be rewritte 
n(t) —— n(t-1) 
tho a 
Or 
MC) eet) 
tM (tM 


Taking the limit as t— o~ 


Or 


>) 


ee ee 


Since the I-P matrix has a 


tne result 


Bye es (Se) 


QO must be true, otherwise there 


which does not make 


This means there are two possible 


R (t) 


1. 


~ 
DL 





as assuming that for "large" t, 
Using the above assumption, 
nas 
E (t-1) 4 
ce ee (eae 
| tt 
(t=) R (t-1)™ 
e-——— 0 + ie 
t co is Eu ai 


results in 


unique inverse (see Theorem 1), 


~ 
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is obtained. The steady state distribution vector is 


achieved by normalizing 7, which results in 


f £cfI-2} (eqn 2.3) 
Bu fi" 





After steady state has been reached it follows that 
the formula 


n(t) = tun LOn et Large” t, 


can be used to approximate the stocks at time tt. The two 
examples from section Al are worked out below to illustrate 


pamaction 2.3 .« 


Example (4) 


This is a hierarchical, four category systen. 
The inital stcck vector is: 
mo) = (0295. 74,2855 11) S 


The transition matrix is: 


.728 .102 0 0 
- .83 .046 0 
a 0 867 .033 
0 0 90 20/ 


The total recruitment is: 


R(t) = B& + (t-1)M where R = 35 and M 


il 
Wn 


The recruitment proportion vector is: 


Cd 


ba (1, Oy 0, 0) e 
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Then 


PPG P2308 0 0 
0 ee -.046 0 
I-P= 
0 0 2.133 =es3 
0 0 0 - 098 
and the inverse is 
Sa 16" 626206 2/63 Zoe 
0 be O62 22034 isis 
(I - p)y-i = 
0 0 72519 Pa el: 
0 0 0 10.204 : 
Using PQUata2en.2.o 
32676 VZ=206 2263 2 OM, 
Sabai (ie, O70. 20) 0 Soo weeeos4 -685 
Ll -°- = eergrFefe,-_o_—_"— 
om 6.902 0 0 Te 519 PF 2s ee 
0 0 0 10.204 


(53, «32, «11, -04) 


is the steady state distributicn vector. 

The ccemputer printout of the terminal session for 
this example, showing the stocks for years 10, 50, 10:03 
900, and 975 can be found in Appendix B, confirming thane 
steady state there will be 553%, 322, 11%, and 4% respec- 
tively, in categories 1,2,3, and 4. As the printout shows, 
Steady state distribution (to two decimal accuracy) 1s 


achieved py t = 975. 


Example (b) 


This 1S a non-hierarchical, three-grade system. 
The inital stcck vector is: 
n(0) = (300, 2007 100) = 


24 


The transition matrix iS: 


380) sal Sy ie: 
P = “/35 «20 0 
-95 0 99 - 


The total recruitment is: 
Rca + (t-1)M where R = 100 and M = 25 . 
The recruitment proportion vector is: 


mea, 007, 223) . 


Then 
fag | =e - <=. 05 
I- P= ~-.75 Aree 0 
-.05 0 - 10 
and the inverse is 
4.267 ao 24133 
(I - Pjp-i = 4 2 2 
Zen 33 4 11.060 : 
Using Equation z.3 
i. 261 -8 Zoi 33 


? (een OT, <23) 
as ee 2 =i) 


\32 
il 
i 
Jes 
ca 
oO 
fh 
r 
ro) 


2.133 .4 11.066 
=e. 43, .09, .48) 


is the steady state distributicn vector. 
The computer printout of the terminal session for 


this example, showing the stocks for years 10, 50, 100, 


350, and 400 can be found in Appendix B, confirming that in 


steady state there will be 43%, 9%, and 48% respectively, in 


categories 1,2, and 3. As the printout shows, steady state 
distribution (to two decimal accuracy) is achieved by 
t = 400 . 


3. Option 3: Multiplicative Recruitment 


This cption addresses the multiplicative increase 
(decrease) in recruitment for each time period. Here, 


Ret-1 , replaces R{t) in Equation 1.2 resulting in 
n(t)> = (t= eee ne at ee (eqn 2.4) 


The values that 6 can assume play an important part in the 
derivation of the steady state equation. First, 98 20 must 
be true, otherwise in even numbered years there would be 
negative recruitment ( R(t)<0 ) which does not make sense in 


manpower models. There are three cases to be examined. 
Case 1 3 Ost o3<1 
Then 
me R(t) = 0 F 
because @f-1—»0Q when O< € <1 Ge Therefore, letting 
t—--~ in Equation 1.2 results in 
isle | 
Oe 
n(I - P) = 0-. 


Since I-P has a unique inverse matrix (see Theorem 1), the 
only solution to the above equation is n= Q. This means 


that suck a system becomes empty in steady state. 
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iim tek (ty) = Roy 


moeh 2S identical to option 1. 


Case 3: e > 1 
Assume that, 





a n (t) a 
limit =— = i 
te o0 S 
exists. TiesSeawounts tO assuming that for “large"™ t, a(t) 
behaves like ett. Using the above assumption, Equation 
2-4 can be rewritten as 
n {t) mele 1) 2 Ret-1 
—_——— =) = —_————..e 1. ae ae ‘ 
ae et-1 ¢@ e* = 


Taking the limit as t—o results in 


Or 
T (SI-P] =Rr . 


Finally, knowing that the matrix @I-P has a unique inverse, 


Since 86 > 1 (see Theorem 1), the result 
@ = Rr (@I-P}-1 


. = ° lay - P : . 
1s obtained. When n 1s normalized, the resulting vector is 


n, the steady state distribution vector where, 


Kkr[{ O1-P j-?t (eqn \2 5) 


‘ a1! 


After steady state has been reached it follows that 


the formula 
n(t) = oD fOr sibarge”  c, 


can be used to approximate the stocks at time t. Two exXam- 


ples are worked out below to illustrate Equation 2.5 


Example (a) 
This is a hierarchical, four category system. 


The inital stcck vector is: 
nq(O) = (129, 74, 28, 11) . 


The transition matrix is: 


ZOO 0 0 
7 0 83 - 946 0 
: 0 0 -86/7 .033 
0 0 0 SS ° 


The total recruitment is: 
R(t) = Ret™-1 where R= 35 , 8 = 1.01. 
The recruitment proportion vector is: 
r= (1, 0, 0770) : 


Then 


20 Caen Ur 0 0 

0 ac ~. 046 0 
6i - P= 

0 0 -143 -.033 

0 0 0 - 108 


and the inverse is 


28 


3. 546 





2.009 - 646 -198 
0 5. 556—m 1-787 ASyite 
[ect — PP j-t= 
0 0 6.993 Die aie 
0 0 0 9.259 
Using Eguation 2.5 
3-546 2.009 ~ 646 -198 
~ 0 5.556 1.787 ~ 546 
a ‘oh (So y(ieeus--0, 0) 
= H = 0 0 6.993 2.137 
oa ni* 2232965 
=a 0 0 0 9.259 
Be 554, <3 04, 


Ss = alee ap 


1s the steady state distribution vector. 
The computer priatout of 
this example, 


the terminal 
Showing the stocks for years 19, 
can be found in Appendix C, 


session for 
there will be 55%, 


50 and 100, 
confirming that in steady state 
31%, 10%, and 3% respectively, in catego- 
maes 1,2,3, and 4, AS tne printout 
distribution (to two 
t = 50 


Shows, steady state 
decimal accuracy) 1S achieved by 
Example (b) 


This is a non-hierarchical, 


three-grade systen. 
The inital stock vector is: 


n(0) = (300, 200, 100) : 
The transition matrix 1S: 
BOO “sao .05 
p= mS * a2 0 0 
05 0 -90 


29 


The total recriudtmeneers: 


R (t) 


= Ret-1 where R 


100 
The recruitment proportion vector is: 


rv=.6 (5 10,. sOnemence nee 
Then 
~ 41 —-.15 =-.05 
@I - P= -.75 at 0 
=205 Q 11 
and the inverse is 
4.03 Oye es. 132 
[Ot =@piet-= 132 oe (coe 1.69 
P53 Poy 9.92 é 
Using Eguation 2.5 
EOS! sp 5 1283 
~ ot (110.0) (270.0207, .23) | 3a73 le oe tee 
l= Sa SS 
= ale 792.15 1.83 ~34 9.92 
= eG 0) cael ag 


The computer 


is the steady state distribution vector. 
this example, 


printout of 


1,2, 


can be found in Appendix C, 
and 36 


confirming that 
9% » 


and 47% respectively, 
As the printout shows, 


the terminal 
Showing the stocks for years 10, 
there will be Wu, 


tion (to two decimal accuracy) iS achieved by t 


session for 


50 and 100, 


in steady state 
in categories 


Steady state distribu- 


30 


109 


= 


Be TOTAL SYSTEH SIZE OPTIONS 


The options in this section are variations of Equation 
1.6 where M(t), representing the change in total systen 
size, is the term that will take on different forms. As a 
result, the steady state equations for these options are 
also obtained in somewhat different forms. In the next two 
subsections the steady state equations will be derived for 


the system size options. 


This option addresses’ the additive increase 
(decrea Se) in total system size for each time period. In 
this case, M(t) = M, a fixed amount of increase (decrease) 


per time pericd. Rewriting Equation 1.6 results in 
n (t) SG OME (eqn 2.5) 


There are 3 cases to be considered. 


Case 1: M< Q 


In this case jobs are being eliminated by a fixed 


amount, M, each year which implies that eventually N(t) = 9 
1s reached in a finite amount of time. This implies that 
n=0i1sS the steady state stock vector. This teans that 


such a system becomes empty in steady state. 


Case 2: M = 0 
When M(t) = O, WN (t) is fixed and can be denoted N. 
It follows that Equation 1.6 can be written as 


n(t) =n(t-1)Q . (eqn 257) 


Since the 


exists {Ref. 4: pp. 40-83, 48-50 Jj, Equation 2.7 can be 


written as 


OE 
n (sO) hee 


Because the first column of I-Q 1S a linear combination of 
all the other cciumas the inverse of I-Q does not exist. 
So, to solve for hn, it is necessary to introduce the addi- 
tional constraint of ni® = N [Ref. 53; p. 30 jz Replacing 
the first column of the 1-Q matrix with a column of ones, 
denoted by (1-Q)*, and replacing the first element of 0 with 


N, denoted 04, tesults tn the following equacion 
ae. (I-Q) * = 0* 


Taking the inverse results in 


Qh = 080 (eos ' (eqn 2.8) 
the steady state stock vector. After steady state is 
reached, then n(t) =n for “large” t. Two examples are 


worked cut pelow to illustrate Equation 2.8 . 
Example (a) 
This 1S a hierarchical, four category systen. 
Tue Aniital Stcex weerome sz 

nO). = (129. ty ea a) : 


The transition matrix is: 


27 ZOCpe el) 2 eno 0 
oe 0 Bis) 5) - 046 0 
0 0 -867 .933 
0 0 0 S02 ’ 
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and 


-898 .102 0 0 
cdo. Os -9 46 0 
ee a 0 1 0) mers 3 
- 098 0 0 = 902 - 


te O* vector 1S; 
O* = ( 242, 0, 9, 0) - 
The recruitment proportion vector is: 
ime O54 Oy 0). « 


Then 


5 WO a 0 
= nl) a ly == 0 46 
ae GO = 168: 
= STs 0 0 


Replacing the first column with ones gives 


1 -~.102 0 
1 ame he) -. 046 
1 - x = 
2 1 0 Sos 
1 0 0 
and the inverse is 
2533 0 lal 


-4.582 3.133 1.084 
meal ele = 6.408 
poo sse20l —1.128 


eo) a. = 


USing Equation 2.8 


om) 
Wel ses, 
LP avs, 
9.824 


uSNe) 0320 2111 S037 
= oG 2 3.133 1.084 365 
=3.394 Yasee 6.408 2.158 
=~50435 ~“=32Z26ulee- | szo 9.824 


[os 
tl 
ame, 
N 
J tay 
A) 
. 
oO 
a 
© 
. 
am 


(128.99, 7is44,. 26. oO pease 
After rounding to the nearest integer the S.S.5.V. 1s: 
n=. (42 9,52 47992 eee 


The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50 and 100, 
can be found in Appendix D, confirming the steady state 
stock vector above. As the printout shows, Steady state is 


achieved by t = 50 . 


Example (4) 


This 1S a non-hierarchical, three-grade systen. 
The inital stcck vector is: 

n(0) = (300, 200, 100) . 
The transition matrix is: 


200 72 ome 


P = Soe DA ve 210 0 
05 ) 29.0 : 
and 
a7 - 164 -096 
Q = Pa se, 2055. 220115 


20S = 0035 ere a> - 


The OS vector 1S: 
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The recruitment proportion vector is: 
Ramet, <0) » soles 


Then 


- 26 en OCs 1G 
fo 70 -= =) he, So eee hc 
=U GD 9 me UF 2089 : 


Replacing the first column with ones gives 


1 -.164 ~-.056 
(I - Q)* = 1 77 Sa0ry 
1 -.004 .0839 


and the tnverse is 


- 430 muo8 2479 
CU = 20) +) aaa ae OC UN ede) 55 16 
= eC Cem Se De 0 69 - 


Using Equation 2.8 


430 090 479 
f=) ( 600, 0, 0) -.611 joie aan 
“4.888  -.981 5.869 


ez o650, 54.4, 287.48) = 
After rounding to the nearest integer the S.5.5.V. 15s: 


mee (258, 54, 287) . 
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The computer printout of the terminal session for 
this example, showing the stccks for years 10, 50 and 100, 
can be found in Appendix D, confirming tne steady state 
stock vector above. As the printout shows, steady state is 


achieved by t = 50. 
Case 3: M > 0 
Equation 1.6 can then be rewritten as 


n(t) a(t) , 2 


if 
N (t) N(t) Ney = 2: 








Iterating the eguation, N(t) = N(t-1) + M, for t = 1,2,.-. 
will shew that, WN(t) = N(O) + tM for all t = lyp2yeee and 


therefore, 








ee M a: N(t- 1) 
lini = =<0 and lini te 
t—-o N(t) t—o N(t-1) + & 
Then wereiteng 
n(t) n(t-1) N(t-1) M 
$$$ = ————__ oe —_____________ 0) + LC 
N (t) N(t-1) [N(t-1) + Gi] Nit)” 
and taking the limit as t—»o yields 
TR =i-0 7 ona Ce 
where 
~ ses n (t) 
n = linit — 
t—- co N (t) 
is the steaay state distribution vector. 
This is identical to the M = 0 case, except here the 
constraint ni* = 1 must hold. Therefore, by replacing the 


first column o£f the I=-Q matrix with a column of Ones 
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denoted by (I-Q)*, and replacing the first element of 0 with 


meamaecnoted 0**, the solution is 


rw 


m= OF#F( (1-0) *]}-? , (om 25 5)) 


the steady state distribution vector. 
ALtecresteadqy state has been achieved, na(t) can be 


found using 
n(t) = OLN(0) + tujr 


fer “large" t. Two examples are worked out below to illus- 


trate Eguation 2.9 . 


Example (a) 
This is a hierarchical, four category systen. 
The inital stcck vector is: 

mcO) =8 (129, 74,828, 11) : 


The transition matrix is: 


£70 Won ie 0 
: 0 .83 .046 O 
0 Bap (ese: 
0 0 0 0D 
and 
-898 .102 0 0 
HiDo 2 OG CO 
as 4 0 NO6Gin 6033 
.098 0 0 apy fF 


wiewuss VeECtOr 1S; 


Oc ea 15-95 07 © 1 e 


at 


The recruitment proportion vector 1S: 


r= (th, 0, Gye) g 


Then 


102 -.102 0 0 
Sipie 217 9=.01c=——n 

a a 0 133s 
-.098 0 0 098 Ya 


Replacing the first column with ones, gives 


1 -.102 0 0 
F i 1 17 * =-.046 0 
nS 1 0 2133 =e 
1 0 0) -098 
and the inverse is 
eo350 320 111 2037 
-4.582 3.133 1.084 365 


I - Q)*)71 = 
=5259 94, 327212 6.408 221455 


=55435 =3e26 | eee 9. 328 : 
Using Equation 2.9 


Ee) Bie. » 320 Slitet tee 
-4.582 3m 138 1.084 5 15S 
=Is354  =BeenZ 6-408 Ze 110 
= 56430 wae Oe ome 9.824 


les? 
iI 
eee, 
—* 
|] 

© 
|] 

© 
|) 

a | 


(5535 23:27 as 


is the steady state distribution vector. 
The ccmputer printout of the terminal session for 
this example, showing the stocks for years 10, 50, 1902 


900, and 950 can be found in Appendix D, confirming that in 
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steady state there will be 53%, 32%, 11%, and 4% respec- 
tively, in categories 1,2,3, and 4. As the printout shows, 
steady state distribution (to two decimal accuracy) is 
achieved by t = 950 . 


Example (b) 
This 1S a non-hierarchical, three-grade systen. 


The inital stcck vector is: 
n(0) = (300, 200, 100) = 
The transition matrix is: 


somes 1D Cli 5 


p= « (5a 20 0 
OS G -90 ; 
and 
- 74 - 164 -096 
Q= + 185 Peo ee C115 


- 085 20055) 2.9115 : 


The 0** vector is: 


The recruitment proportion vector is: 


r= (270, .07, .23) = 


Then 
a2 0 emo dae. 96 
Os 5 ahs 2797 2 =.012 
Seieisy ey OO) 2059 : 


Replacing the rirst column with ones, gives 


~~ 


1 -.164 -.096 
(I - Q)* = 1 mo =2012 
17 -.004 #£.089 


a0 


and the inverse is 


Po CaO 2 8 79 
((2 = oy )=) = [ete is tices ic 
STi eee = ches. See 2 


Using Eguation 2.9 


430 .091 -479 
T= (17°07 -.611 1.127 -.516 
-4.888 -.981 5.869 


= (5.43. <O975. 45) 
is the steady state distribution vector. 
The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50, 100, 
350, and 400 can be found in Appendix D, confirming that in 


steady state there will be 43%, 9%, and 48% respectively, in 


categories 1,2, and 3. As the printout shows, steady state 
distribution (to two decimal accuracy) is achieved by 
t= 400 . 


This section addresses the option of multiplicative 
increase idecrease) in total system size for eaca tine 


Feriod. Starting with 


N (t) ON (t-1).,. “sor £22 
it follcws that 

N(t) = eN(0), fort = 1,2,... 
from which M(t) can be determined as 


~ 


M(t) = N(t) - N(t-1) = (6-1) e6%-1N(0) ~ 
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Then, letting (6-1)e@%-!1N (0) replace M(t) in Equation 1.6, 


results in 
mat) = n{t-1)0 + [ (@-1)e°-IN()j] rc - (eqn 2. 10) 
First, ©6> 0 must be true otherwise negative or zero system 
size may result. This leaves two cases to be examined. 
Case 1: 0O< €é€< 1 
In this case 
limit N(t) = 0 
t ——~ 0 
because aS t—e oo , Ot-!—~» 0 for 0 <e<1. Finally, if 
8 = 1 then 6-1 = 0 and N(t) = 0 for all t. Therefore, 
nh = limit n(t) =9 
t > © 
must also hold. This means that such a system becomes empty 
in steady state. 
Case 2: e@> 1 
Equation 2.10 can be rewritten as 


n(t) | p(t-1) oN) 


| 
N (t) N{t) N(t) = : 


Since N(t) can kre written as ON(t-1) or as e°N(9), the above 


equation gives 





n (t) n(t-1) 1 (9-1) o*-1N (0) 
eas eee ee + SS a © 
N (t) N(t-1) 0 6*N (0) = 
Or 
n (t) n(t-1) 1 e-1 
i ee ge eQ + Cc = 
N ({t) N(t-1) 6 Se oad 


Taking the iimit as t— » eo results in 
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les 2 
iN 
Iss 





or 


[2 


[Gl=O]% = Cease 


where 
n (t) 
t— oo N (t) 





is the steady state distributicn vector. 
Since 6 > 1 in this case, the matrix 98I-Q has a 


unigue inverse according to Theorem 2, so the result 


Cand 


N= (S51) Sr es Seria (eqn 2.11) 


is obtained for the steady state distribution vector. After 


steady state has been reached it follows that the formula 


n(t) = @°N(0) @ 


can ke used te find the stocks at time t. Two exampies are 


worked cut below to illustrate Equation 2.11. 
Example (a) 
This 1S a hierarchical, four category system. 
The inital stcck vector is: 

n(0) = (429.. 749 26,8?) : 


The transition matrix is: 


i 0 -83 .046 9 
7 0 0 Ney 4088 
0 0 0 Oe 
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and 


-898 .102 0 0 
SPAS eS: - 0 46 0 
: al 0 -867 .033 
- 098 0 0 702 - 


The multiplicative factor is: 
Guten 30 1 y< 


The recruitment proportion vector is: 


Then 
~ 112 -. 102 Q ¢) 
-. 124 ~-18 -.046 Q 
6I - Q = 
-. 1 ¢) ma 3) =2 033 
-. 098 0 0 - 108 


and the inverse is 


Soe S440 M10. 1G 3. 09 
51.04 34.48 11.09 eee 
SiG 3)5) Atcha SS) a Were an) 4.94 
50.28 28.49 Beli 1 Zeruo . 


Hoga “O)ig*- = 


Dsing Bgquation 2.11 


552.41 31.40 10.10 3.09 
51.04 34.48 11.09 SERS 
DU Se ice — 16. 17 4.94 
02 Zo Zo. 49 Peivee ic. OG 


[2 


(-01) (1,0,0,0) 


HI 


(6554, .314, .101, ~031) : 


is the steady state distribution vector. 
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The computer printout of the terminal session for 
this example, Showing the stocks for years 10, 50, and 100, 
can be found in Appeniix E, confirming that in steady state 


there will be 55%, 31%, 10%, and 3% respectively, in catego- 


ries 1, 2,3, and 4. As the printout Shows, steady state 
distribution (to two decimal accuracy) 1S achieved by 
t = 50 =. 


Example (b) 


This iS a non-hierarchical, three-grade systen. 
The inital stcck vector is: 

n(0) = (300, 200, 100) . 
The transition matrix is: 


= Om ts ec SD 


p= et Se 0 0) 
0S Q 90 P 
and 
: a 78 164 996 
9 = . OD a2 03 See lis 


s 005 ~ 0035 wee 145 : 
The multiplicative factor is: 
G = 1.01 . 


The recruitment proportion vector is: 


Then 


-27  -.164 -.096 
-.785  .8065 -.0115 
-.085 -.0035 .29985 


Gl = Ces 
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and the inverse is 


eS 2 9.41 45.27 
44.69 10.52 44.79 
40.70 8.50 50.81 - 


(@I - Q)-? = 


Using Equation z.11 


45232 92-41 45.27 
44.69 10.52 44.79 
40.70 8.50 50.81 ° 


(a0 4) (= 7, Ose. 23) 


[32 
T 


(-44, -09, 47) 


is the steady state distribution vector. 

The computer printout of the terminal session for 
this example, showing the stocks for years 10, 50, and 100, 
can be found in Appendix E, confirming that in steady state 
there will pe 44%, 9%, and 47% respectively, ir categories 
wee, and 3. As the printout shows, steady state distribu- 


tion (tc two decimal accuracy) is achieved by t = 100. 
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Iii. SUMMARY 
This chapter is a synopsis of the options that describe 
the five different conditions of growth. In all option? 
the restriction of w, > 0 for alli = 1,2,...,k was assumed, 
Since attrition must be allowed from any category in a 
Manpower systen. Options 1 through 3 employed the kLasic 
prediction equation (Eguation 1.2) and used constant, addi- 


tively and nultiplicatively increasing (decreasing) recruit- 


ment, R(t), to derive the steady state behavior of a 
manpower system. For options 4 and 5, Equation 1.6 used 
additively and muiltiplicatively increasing (decreasing) 


total system size, M(t), to derive the steady state behavior 
of a manpower system. 

Derivation of the steady state stock vector or distribu- 
tion vector required some standard results of matrix algebra 
Summarized in Theorems 1 and 2 of Chapter II. In addaiwre 
to the analytical derivation of the steady state equations, 
examples have been worked out in each case (excepting some 
trivial cases) to illustrate the use of the analytical 
results. The formulas have been programmed in the APL 
language and the program listing is given in Appendix F. 
This program (called a function in APL) has been integrated 
into the main APL program at NPS which forecasts future 
Stocks of a Manpower system using the Markov Chain theory 
Summarized in Chapter tI. Appendices A through E£ Show the 
printouts of the terminal sessions where the examples 
referred to akove are worked out to verify the steady state 


results obtained in each analytical example. 


w 
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APPENDIX A 


COMPUTER PRINTOUTS OF EXAMPLES FOR OPTION 1 


START 
DO YOU WISH TC ENTER DATA? 
oo NC 
1 YES 
: 1 
ENTER THE NUMEER OF THE MODEL TYPE 
1 MARKOV HIERARCHICAL 
2  MAKKOV LENGTH OF SERVICE 
3. MARKOV GENERAL 
4 VACANCY 


129 74 28 11 

ENTER THE PROMOTION RATE VECTOR. 
THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE 

3 CLASSES. 

: MiG2eeOuG 2033 
ENTER THE WASTAGE RATE VECTOR. 
THE 4TH VALUE SHCULD INCLUDE THE TOTAL RATE FOR THAT 
STATE DUE TO EITHER PROMOTION OR WASTAGE. 

; migee tu 27.098 
ENTER THE NUMBER OF THE RECRUIT TYPE 

1 FIXED RECKUIT VECTOR | 

2 ACDITIVE (RECRUIT S E). 

3 MULTIPLICATIVE (RECROL SIZE) 

ri ADDITIVE (SYSTE E) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

. 1 
ENTER R (RECRUITMENT VECTOR) 

; 3500 0 
ENTER THE PERCENT CODE 

6) NO GRADE PERCENTAGES 

1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZsS 

2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 

q OUIT PROGRAM 


1 
ea ac Pre elOmseh Lib ENTERED DATA? 
1 LES 
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e tO 
CO =i 
WOrsd 

ho 


t3— Cory : 
© 
wad $2) 


ON tlonroooow7 
© 
8) 


io Ran LIKE TO CHANGE ANY OF THE DATA? 


1 YOTS) 
i QUIT PRCG KAM 


0 
ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
: 10 
DO YOU WISH TO SEE THE INTERVENING YEARS? 
N 


0 

i N PERCENT R 
0 1 129 33) 
2 74 31 
3 28 ir 
U4 11 5 
TOTAL 242 100) 
10 1 129 33) 
2 ng 32 
3 27 17 
4 10 4 

TOTAL DAD 100 35 


— Se eo 


OT ae ee TO SEE ANY OTHER YEARS? 
N 


71 WES 
7 QUIT PRCGRAM 


1 
ENTER THE NUMEER OF THE YEAR YOUR Wise em oe 
50 
DO YOU WISH TC SEE THE INTERVENING Yea Rog 
1 YES 


0 
50 1 129 53 
2 Da 33} 
3 27 n 
Y 9 4). 
TOTAL 242 100) 35 
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o, PCS SEE ANY OTHER YEARS? 


1 ek 
U QUIT PROGRAM 


1 
ENTER THE NUMBER OF THE YEAR YOU WISH 10 SEE 
: 100 
DO YOU WISH TC SEE THE INTERVENING YEARS? 
1 Yes 


0 
100 1 129 33| 
2 Oy a2 
3 Fou 11 
+ 3 i 
TOTAL 242 100 35 


DO FO ey al TO SEE ANY OTHER YEARS? 


1 EES 
7 QUIT PRCGFRAM 


. 0 
DO, YOU WISH TO SEE THE STEADY STATE VECTOR? 
am ves 


999 1 129 53 
2 77 32 
5 27 11 
rm 9 4 
TOTAL 22 100) 35 
ARE YOU THROUGH? 
0 NC 
1 YES 
' 1 
2- Example (b) 
START 
DO YOU WISH TC ENTER DATA? 
0 NC 
1 YES 
: 1 
ENTER THE NUMBER OF THE MODEL TYPE 
1 MARKOV HIERARCHICAL 
2 MARKOV LENGTH OF SERVICE 
3 MAKKOV GENERAL 
Uy VACANCY 
. 3 
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ENTER N(ENT TIAL STOCK VECTOR, 


300 200 100 
foe Ce MATRIX) BY ROWS 
ER 11TH ROW 


E RECRUIT TYPE 
CTO 

T S E) 

RECROL SIZE) 


YSTEM 
(SYSTEM SIZE) 


ee 


H 
E 
au 
41 


F 
ENTER R (RECROITMENT VECTOR) 


EO baa 
ENTER THE PERCENT CODE 
QO NO GRADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GKADE SIZE 
2 GRADE SIZE AS PEPCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 


NO 
1 YES 


Ps 


ee 
nit 


7) 
heJededtj fy—ad 

Oo 

Oro 


Oil 
O 
6) 


U LIKE TO CHANGE ANY CF THE DATA? 
S 
I 


ormogd 


T PRCGRAM 


CO OKZK WH OmoOooeor -— 


ENTER THE NUMBER OF THE YEAR YOU WISH 10 SEE 
: 10 

DC YOU WISH TC SEE THE INTERVENING YEARS? 
1 ¥Es : 
i 0 


50 


ie N PRHGENT R 


0 1 300 50) 
Z 200 3a 
5 100 17 
TOTAL 600 100 
10 1 351 48 
2 75 10 
3 3).9) 2 41 
TOTAL 728 121 1090 


DO ee ZC SEE ANY OTHER YEARS? 
Oars) 
QUIT PRCGRAM 


~] at © 


ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
j 50 
DO YOU WISH TO SEE THE INTERVENING YEARS? 

1 YES 


0 
50. 


100 


oe a TC SEE ANY OTHER YEAKS? 


1 OLS 
7 QUIT PRCGRAM 


1 
ENTER THE NUMEER OF THE YEAR YOU WISH T1O SEE 
1 00 
DO YOU WISH TC SEE THE INTERVENING YEARS? 
1 YES 


0 
100 6 4 3 
9 S 
8 48 
3 146 100 


a ol ies oun ANY SOLER eA RS © 


1 ier 
i) QUIT PRCGRAM 


0 
ooo oN 1G OBE TAS oP eAD ot ATE VECTOR? 
1 toro 
1 
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100 


| 


1 
2 
3 


TOTAL 


a9 


ARE YOU THROUCH? 


a2 


COMPUTER PRINTOUTS OF EXAMPLES FOR OPTION 2 


1. Example (a) 


START 
DO YOU WISH TO ENTER DATA? 
0 NC 
1 Vers 
1 
feok THE NOMEBER OF THE MODEL TYPE 
1 MARKOV HIERARCHICAL 


2 MARKOV LENGTH OF SERVICE 
m MARKOV GENERAL 


VACANCY 
; 1 
ENTER N(INITIAL STOCK VECTOR) 
: 129 74 28 11 
ENTER THE PROMOTION RATE VECTOR. 
THIS VECTOR SHOULD INCLUDE THE PROMOTICN RATES FOR THE FIRST 
3 CLASSES. 
: P10 we 046. 2033 
ENTER THE WASTAGE RATE VECTOR. 
THE 4TH VALUE SHCULD INCLUDE THE TOTAL RATE FOR THAT 
STATE DUE TO EFITHER PROMOTION OR WASTAGE. 
: -17 .124 .1 .098 
ENTER THE NUMBER OF THE RECRUIT TYPE 
1 FIXED RECRUIT VECTOR 
2 ADDITIVE (RECKULT SIZE) 
3 MULTIPLICATIVE {kECRUI SIZE) 
4 ALDITIVE (SYSTEM SIZE) 
5 MULTIPLICATIVE (SYSTEM SIZE) 
: 2 
R (RECRUITMENT VECTOR) 


ENTER 


5 
ENTER THE PERCENT CODE 
QO NO GRADE PERCENTAGES 
1 GPADE SIZE AS PEECENT OF TOTAL GRADE SIZE 
2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM . 
1 
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Sonar YOU LIKE TO SEE THE ENTERED DATA? 


TOTAL RECRUITMENT 
jRECRUITWENT PROPORTION VECTOR 


Mm won OOOoe 


Ww HO 


aes Ne LIKE 10 CHANGEV ANY) OF ei iih eos a 


1 GES 
7 QUIT PRCGRAM 


0 
ENTER THE NUMBER OF THE YEAR YOU WISH 10 SEE 
. 10 
DO YOU WISH TO SEE THE INTERVENING YEARS? 
1 Ys 


0 
T N PERGENG R 

0 1 129 a8. 
2 74 31 
3 28 12 
4 11 5 
TOTAL 242 100) 
10 l 248 $3} 
2 106 Aad 
3 29 7 
4 19 3) 

TOTAL 392 162 80 


Bie cues rh ca TO SEE ANY OTHER YER S. 
: OuTT PRCG RAM 
1 
ENTER THE NUMEER OF THE YEAR YOU Wish ei08s57. 
30 - 
DO YOU WISH 710 SEE THE IN TERVENTIG Ete 
1 Yes 
0 


54 


50 1 980 98 
2 325 an 
3 jez S 
4 39 z 
TOTAL too’ 700 280 


aoe NG Pee oth ANY SOLa en LEAR SZ 


1 HOES 
7 QUIT PRCGKAM 


; 
ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
; 100 
OS dat TO SEE THE INTERVENING YEARS? 
1 YES 


° 
@ 


0 
100 


} 530 


= eS Se ee cee a eS eS ee eee 


Oo i TC. SEE ANY OG2HER YEARS: 


1 YES 
if QUIT PRCGEAM 


4 
feean CHE NUMEER OF THE YEAR YOU WISH 0 SEs 
9 00 
eeyo WiSH TO SEE THE INTERVENING YEARS? 
1 Yes 


0 
900 


WN = 


TOTAL 3 


) 4530 


no seo TC SEE ANY OTHER VERB SZ 


1 YES 
7 QUIT PRCGKAM 


4 
Peon tHE NUMBER OF THE YEAR YOU WISH TO SEE 
S) 18. ” 
POmroumnhEoH IC SEE THE INTERVENING YEARS? 
1 YES 
0 
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2). 
32 
11 
i} 

TOTAL 3 13830) 4905 


DOR aoe TC SEE ANY OTHER YEARS? 


1 LES 
7 QUIT PRCGRAM 


0 
DO YOU WISH TC SEE THE STEADY (5 Ware clon. 
NC 
1 TES 


4905 


ARE YOU THROUGH? 
0 NO 
1 YES 
1 
2- Example (b) 


START 
DO YOU WISH TC ENTER DATA? 
0 NO 
1 YES 
: 1 
ENTER THE NUMBER OF THE MODEL TYPE 
1 AARKOV HIERARCHICAL 
2 MARKOV LENGTH OF SERVICE 
3 MARKOV GENFRAL 
ri VACANCY 
: 3 


ENTEP N (INITIAL STOCK VECTOR) 
300 200 100 
FNTER P (TRANSITION MATRIX) BY ROWS 
ENTER 17H ROW 
: ~6 .15 .05 
ENTER 2TH ROW : 
: a (is eee 
ENTER 3TH ROW 


sey SO mee 
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ENTER THE NUMBER OF THE RECRUIT TYPE 
IXED RECRUIT VECTOR 
2 ADDITIVE (RECRUIT SIZE) 
3 MULTIPLICATIVE (RECRUL SIZE) 
4 ADDITIVE (SYSTE SIZE) 
5 MULTIPLICATIVE (SYSTEM SIZE) 
: 2 
ENTER R (RECRUITMENT VECTOR) 
° 707 23 
ENTER ADDITIVE INCREASE 
: 25 
ENTER THE PERCENT CODE 
O NO GRADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GRADE SI 
2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE 
7 QUIT PROGRAM 
; 1 
HOULD ¥ Ou LIKE TO SEE THE ENTERED DATA? 
1 YES 
; 1 
P MATRIX 
0.6 0.15 0.05 
M275 0.2 0 
0.05 0 0.9 
| N VECTOR 
300 200 100 
OPTION =2 
INC =25 
TOTAL RECRULIMENT 
100 
RECRUITMENT PROPORTION VECTOR 
pee 0-07 0. 23 
ROULD YOU LIKE TO CHANGE ANY OF THE DATA? 
N 
7 Y ES 
7 QUIT PRCGRAM 
; 0 
ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
. 10 
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DO YOU WISH TC SEE THE INTERV ENG Ea ko 


NO 
YES 
; 0 
TN PERCENT — BL 
0 1 300 HO 
Z 200 33 
3 100 17 
‘ BOTA LT 600 100 
10 1 801 53) 
2 162 11 
16; | 1} 
TOA i 1523 254) 325 


DO etary TC SEE ANY OTHER YEARS? 


7 Ort PRCGERAM 
1 
ENTER THE NUMBER GF DHE YEAR YOU CWEZSH TO "SEE 
; 50 
DO YOU WiSH TO SEE THEVEN TER VENEGAS 
om 


0 
50 6 
Q 
> 
1 18-25 


CS ae ews th ANY OTHER Years. 


1 YES 
7 QOUIT PROGRAM 


1 
ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
: 100 
DO, You WISH TC SEE THE INTERVENING YEARS? 
NO 
1 YES 


0 
100 1 9066 4 4 
2 1903 y, 
3 9469 46 
TOTAL 20438 3406 2575 


OOF ae TC SEE ANY QUHER Yiato ce 


1 iS 
if QUIT PRCGKAM 


1 
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ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
: 350 

DO YOU WISH TO SEE THE INTERVENING YEARS? 
Ys 


0 
350 1 32550 43) 
2 0853 9 
3 35585 47 

TOTAL 74987 12498 8825 


ot DC SEGeaANY OTHER YEAKS<? 


1 YES 
7 QUIT PRCGRAM 


1 
ENTER THE NUMBER OF THE YEAR YOU ¥ISH TO SEE 
; 4 00 
DOLYOU WISH TC SEE THE INTER VENGNG YEARS? 
1 YES 


te 


0 
400 1 37246 4 3 
2 7843 a 
3 40808 4S 
TOTAL 85897 14316 10075 


DO Sain me SEE ANY OTHER WEARS 2 


1 Ob ips: 
7 QUIT PROGRAM 


0 
DO YoU WISH TC SEE THE STEADY STATE VECTOR? 
N 
1 YES 
: 1 


Meme aNTAGES ARE IN STEADY STATE 
poo 1 37246 43 
Z 7843 2 
3 40808 48 
TOTAL ao 897 100) 10075 


ARE YOU THROUGH? 
0 ° Ne 
1 YES 
: 1 


APPENDIX C 
COMPUTER PRINTOUTS OF EXAMPLES FOR OPTION 3 


1. Example (a) 


START 
DO YOU WISH TO ENTER DATA? 
0 NC 
1 YES 


j 


ENTER THE NUMBER OF VITHE Mepel Iyer 
1 MARKOV HIERARCHICAL 
2 MARKOV LENGTH OF SERVICE 
3 MAKKOV GENERAL 
4 VACANCY 


1 
ENTER N (INITIAL STOCK VECTOR) 


129 74 28 11 
oe PROMOTION RATE VECTOR. 
S 


: ae SHOULD INCLUDE THE PROMOTICN KATES FOR THE Shae 


- 10 
ER THE LO 

UTH VA E THE TOTAL RATE FOR THAT 
E DUE ON OR WASTAGE. 


Ss 
Esa) 
mS 
Ot'rs 
raM 
ah, 


ti 
<= 


(n= LINO + 
tz4 
ru 


e 
iS 
1 
Hi 
U 
1 
Hq 
ag 
D 
U 
C 
U 


br <td <jto fN 
be4$ Kh ETI 4 OO 


(RECRULTAENT VECTOR) 
me O250 
ULTIPLICATEVE FAGIOR 


Vigo 
So 


raj rr 
Org 


OTAL GRADE 
RIGINAL GR 


-~ 


ST. 
ADE E 


Mw 


OAD” 
a COMMON = SW WW repre maye 


60 


eae NO Deke TO SEE sEaE ENTERED QaTraA ? 
1 YES 


@ 

i Aa 8) 

Se) 

CO 

> om 

9 
CO = 
WIOrs 

NO 


OHO 
HH = 
Ox] 
=a=e< 
ty 
lh lrmoAGDoovorn 
ml Corl 


TOTAL RECRUITMENT 
ee FROPORTION VECTOR 


Mm BON GOOO 


Ww MO 


WOULD YOU LIKE TO CHANGE ANY CF THE DATA? 

7 ort PROGRAM 
; 0 

ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
10 | 

DO YOU WISH TC SEE THE INTERVENING YEARS? 
1 ¥8s 


0 
T N PERCENT R 

0 1 129 ( 53 

2 74 31 

5 28 12 

m 11 5 

TOTAL 242 100) 

10 1 oa 54) 

2 79 an 

3 2a in 

i 0 i 
TOTAL 253 105 38 


ae ee mO. SEE ANY OTHER YERe se 


1 eo 
| QUIT PRCGRAHM 


1 
Beer THE NUMEER OF THE YEAR YOU WISH TO SEE 
2) 0, 
DO YOU WISH TO SEE THE INTERVENING YEARS? 
a He. 
0 


61 


50 


Fw 


TOTAL 365 1 | S|) 


= Sap cep cae Ge SEP UD Sc cee ce cee OE ee ee ee ee ee eee ee eee ae 


DO, YOU #ISH TC SEE ANY OTHER YEARS? 


1 YanS 
i QUIT PRCGEAM 


ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
° 100 

DO. YOU WISH TC SEE THE INTERVENING YEARS? 
1 YES 


0 
100 1 336 55 
2 190 24 
3 61 10 
4 19 3 
TOTAL 606 250 94 
DO YOU WISH TC SEE ANY OTHER YEARS? 
1 YES 
7 QUIT PRCGRAM 
0 | 
DO YOU HISH TO SEE THE STEADY STATE VECTOR? 
“ | 
1 Y ES 
‘ 1 
PERCENTAGES ARE IN STEADY STATE 
999 1 336 55 
2 190 31 
3 61 10 
4 19 5 
TOTAL 606 109% 94 
ARE YOU THROUGH? 
NO 
1 YES 
: 1 
2. Example (b) 
START : 
DO YOU WISH TO ENTER DATA? 
0 NC 
1 YES 


1 
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ENTER THE NUMEER OF THE MODEL TYPE 
1 MARKOV HIERARCHICAL 
2 MARKOV LENGTH OF SERVICE 
3 MARKOV GENERAL 
4 VACANCY 
: 3 
ENTER N (INITIAL STOCK VECTOR) 
i 300 200 100 
ENTER P (TRANSITION MATRIX) BY ROWS 
ENTER 17H ROW 
; ~6 215 2.05 
ENTER 2TH ROW 
i P75 22.0 
ENTER 3TH ROW 
j -05 0 .9 
ENTER THE NUMBER OF THE RECRUIT TYPE 
1 FIXED RECRUIT VECTOR 
2 ADDITIVE (RECRUIT SIZE) 
3. MULTIPLICATIVE (RECRUL SIZE) 
4 AOL IAEA SY STE SIZE) 
5 MULTIPLICATIVE (SYSTEM SIZE) 
F > 


ENTER R (RECRUITMENT VECTOR) 


707 23 
ENTER MULTIPLICATIVE FACTOR 
; 1.01 
ENTER THE PERCENT CODE 
Q NO GRADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 
2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
i QUIT PROGRAM 


1 
mia ate EIKE TO SEE DRE) ENTERED SAMA < 
1 vos 


POTAL RECRUITMENT 
100 


RHGEULTMENT PROPORTION VECTOR: 
wee. 07 0. 23 
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he pee LIKE TO CHANGE ANY OF THE DATA? 
7 OurT PROGRAM 
0 
ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
fe 
mae) Speke eagua TO SEE THE INTERVENING YEARS? 
1 YES 


Q 
ee 21 
0 1 300 30 
2 200 33 
3 100 17 
TOTAL 600 100 
10 1 370 49 
2 78 19} 
S 302 4 
TOTAL 761 127) 109 


DO + OO aoe fC she ANY OTHera yon 52 
I 


] OaTT PRCG RAM 
1 
ENTER THE NUMEER CF THE YEAR YOUGWISHRIGeSsEe 
5 0 
Te eee TG SEE THE INTERVENING YEARS? 
UC eS 


0 
50 1 575 
2 121 
603 
163 


eC ee DE SEE ANY OTHER YEARS? 


1 Oa 
7 QUIT PROGERAM 


1 
ENTER THE NUMBER OF THE YEAR YOU WIESE 20 See 
100 
DO YOU WISH TO SEE THE INTERVENING YEARS? 
1 Yes 7 
: 0 


64 


268 


= ee ee ee ee ee ee ee ee ee es ee ee eee eo ee eo eee ee ee es eee 


ae a TC SEE ANY OTHER YEARS? 


1 YES 
7 QUIT PRCGRAM 


0 
DO ee er? i DOMSsEEeIiE SLEADY STATE VECTOR? 


FERCENT AGES ARE IN STEADY STATE 
399 1 947 44) 
2 199 3 
3 996 4 7 
TOTAL 2141 100) 268 


NC 
YES 


4 


Sa YOU THROUGH? 


APPENDIX D 
COMPUTER PRINTOUTS OF EXAMPLES FOR OPTION 4 


1. Case 1 Xxample (a) 


START 

DO YOU WISH TC ENTER DATA? 
0 NO 
1 bg as 


OF feb L iYre 
ARCHICAL 
Lh OP @sERVICE 


ENTER N (INITIAL STOCK VECTOR) 
29 74 28 11 

HE PROMOTION 

CTOR SHOULD 
SES. 


RATE VECTOR. 
INCLUDE THE PROMOTION RATES FOR THE Site 


tite SB wm ike 


Pas 
we 
© 
NO 


ATE FOR THAT 


omre 
ad 
Mitt =) 
ad— =Oape 
tr 


reg it < oish 
© 


(nf WN 
re Sate 
CSmWOOKOHms G Tas 
HForPomxim~ ma 
Hr Hejeitd 
Kd ee 
FO HIS PY 
kf ort H4 bei tf tr 

DM Atty 


<= 


ENTER RPROP (RECROULTMENT PROPOKTEON. VEGIOR) 


bees G0 
ENTER ADDITIVE INCREASE 


0 
ENT oR THE PEREENT CODE 
O NO GRADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GEADE Size 
2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 3 
1 


ee “ LIKE TO SEG THE. ENTERED Gee. 
1 YES 
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P MATRIX 
0.728 0.102 0 0 
0 0.83 0.046 O 
0 0 0.867 0.033 
0 0 0 0.902 
N VECTOR 
W574 28 11 
OPTION =4 
INC =0 
RECRUITMENT PROPORTION VECTOR 
mo 0 0 
WOULD ¥ ou LIKE TO CHANGE ANY CF THE DATA? 
1 YES 
7 QUIT PROGRAM 
; 0 
ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
j 10 


DO YOU WISH TC SEE THE INTERVENING YEARS? 
1 YES 


0 
T N PERCENT R 
0 1 129 53 
2 74 37 
3 20 12 
4 11 5 
TOTAL DNS 100) 
10 1 129 53) 
2 Ti 32 
3 2) 11 
4 10 4 
TOTAL 242 100 35 


a 2 a ee ee ee es 


eo me SEE ANY OTHER YEARS? 
] OurT PRCG EAS 
1 
vols NUMBER OF THe YEAR YOU WISH TO SLE 
50 
mGey OU WISH TO SE& THE INTERVENING YEARS? 
1 YES 


She 


DO YOU WISH TO SEE ANY OTHER YEARS? 
0 NC 
1 YES 
7 OUIT PRCGRAM 
1 
ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
: 100 
DO. YOU WISH TO SEE THE INTERVENING YEARS? 
1 YES 


Q 
100 1 1238 33 
2 Od 32 
3 27 11 
4 3 + 
TOTAL 242 100 35 


ee Nicaea [OC SEERA Y*OT HER =MEAnS < 


1 Yes 
7 QUIT PRCGKAM 


e 
@ 


0 
DO YOU WISH TO SEE THE STEADY STATE VECTOR? 
1 YES 
. 1 
PERCENTAGES ARE IN STEADY STATE 
999 1 129 53) 
2 a9, 32) 
3 ph) 11 
4 9 4) 
TOTAL Di 100 35 
ARE YOU THROUGH? 
NO 
1 YES 
‘ 1 
2. Case 1 sxample (b) 
START 
DO YOU WISH TC ENTER DATA? 
0 NC 
1 YES 
: i. 
ENTER THE NUMEER OF THE MODEL TYPE 
1 MARKOV HIERARCHICAL 
D MARKOV LENGTH OF S™®RVICE 
S MARKOV GENERAL 
4 V ACANCY 
; 3 
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ENTER N(INETTIAL STOCK VECTOR) 
300 200 100 


ENTER P (TBANSITION MATRIX) BY ROWS 
ENTER 11H ROW 
j 6 2.15 .05 
ENTER 2TH ROW 
J AG 
ENTER 31H ROW 
: _ Sy tay fe 
ENTER THE NUMPER OF THE RECRUIT TYPE 
1 FIXED RECRUIT VECTOR 
2 ADDITIVE (RECRUIT SIZE) 
3 MULTIPLICATIVE {RZCRUI SIZE) 
4 ADDITIVE (SYSTEM SIZE) 
5 MULTIPLICATIVE (SYSTEM SIZE) 
' mM 
ENTER REROP (RECRUITMENT PROPORTION VECTOR) 
' a Og eee 
FNTER ADDITIVE INCREASE 
: 0 
ENTER THE PERCENT CODE 
0 NO GRADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 
2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGR 
‘ 1 
WOULD ¥ ou LIKE TO SEE THE ENTERED DATA? 
1 YES 
: 1 
P MATRIX 
0.6 0215 0.05 
O75 0.2 0 
0.05 0 Care 
N VECTOR 
300 200 100 
CETION =4 
INC =0 
RECRUITMENT PROPORTICN VECTOR 
fey 0.07 0.23 
ae Noe LIKE TO CHANGE ANY OF THE DATA? 
1 YES 
7 OuULIT PRCGRAM 
e 0 ; 
ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
° 10 
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CO TC SEE THE INTERVENING YEARS? 
Ae aes 


0 
_ 7 NN ER eee 
0 1 300 S00 
2 200 | i 
3 100 17 
TOTAL 600 100 
10 1 284 47 
2 61 | i 
3 OAS 4 3 
TOTAL 600 100 We 


ed eee TO SEE ANY OCTEER yeas 


1 NGES 
7 QUIT PRCGRAM 


4 
ENTER THE NUMEER OF THE YEAR YOU MESH 10552. 
ay) 
DO YOU WISH TO SEE THE INTERVENING YEARS? 
1 YES 


0 
50 1 258 43 
2 54 2 
3 287 054 
TOTAL 600 100 69 


De TC SEE ANY OTHER YEARSS 


1 Ls 
i GUIT PRCGRAM 


4 
BNTER THE NUMEER CF THE YEAR YOUU sheer oc. 
100 
was ae - eanl TC SEE THE INTER VENENGe ye Ao 
1 Des: 


0 
100 1 258 43 
2 54 3 | 
3 287 48 
TOTAL 600 100). 69 
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a ee POSSE ANOTHER YEARS? 


1 YES 
i} QUIT PRCGRAM 


0 
DO YOU WISH TO SEE THE STEADY STATE VECTOR? 
1 YES 
1 
PERCENTAGES ARE IN STEADY STATE 
999 1 258 43 
S4 9 
3 287 48 
TOTAL 600 100 69 
ARE YOU THROUGH? 
0 NC 
YES 
; 1 
3. Case Zz Example (a) 
START 
DO, You WISH TC ENTER DATA? 
1 YES 
ENTER THE NUMBER OF THE MODEL TYPE 
1 MARKOV HIERARCHICAL 
2 MARKOV LENGTH OF SERVICE 
5 MARKOV GENERAL 
4 VACANCY 


: 
ENTER N (INITIAL STOCK VECTOR) 


129 74 28 11 


ENTER THE PROMOTION RATE VECTOR. 
THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE 

3 CLASSES. 

: -102 .046 .033 
FNTER THE WASTAGE RATE VECTOR. 
THE OTH VALUE SHOULD INCLUDE THE TOTAL RATE FOR THAT 
STATE DUE TO EITHER PROWOTION OR WASTAGE. 

‘ -17 .124 .1 .098 
ENTER THE NUMBER OF THE RECRUIT TYPE 

FIXED RECRUIT VECTOR z 

DePASDITIVE (RECRUID STZE) 

3 MULTIPLICATIVE (RECRUIT SIZE) 

4 ADDITIVE (SYSTEM SIZE) 

5 MULTIPLICATIVE (SYSTEM SIZE) 

; 4 
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Ja TE se) IE 


ENTER RPROP (RECRUITMENT PROPORTION VECTOR) 


E T2000. 6 
ENTER ALDDITIVE INCREASE 
5 


ENTER THE PERCENT CODE 
NO ADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 
2 GEADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 


1 YES 


UITMENT PROPORTION VECTOR 

U LIKE TC CHANGE ANY OF THE DATA? 

S 

IT PRCGRAM 

ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
i 10 


ae hee ae TO SEE THE INTERVEM Geyer. 
1 ees 


v 
T N PERC rh R 
0 1 129 J 
2 74 a4 
3 Ze 12 
+ 11 3 
TOTAL 242 100) _ 
10 1 165 364 
Ze 89 30 
5 28 10 
4 10 3) 
TOTAL Z92 121 4 7 


OS ee POR IEee ANY G@THER YEARS? 
‘ 0 Ore PRCG RAM 
ENTER THE NUTPBER VOR THES YEAR YOU WESH TO SEE 
5 0 
PeerCU WISH TO SEE THE INTERVENING YEARS? 
1 YES 


0 
50 


WN) 


596 
24 
i 
TOTAL 492 208 


a a ee ee ee ee ee co 
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an ee LCe ose COOLER YEAKS< 


1 YES 
7 QUIT PKCGRAM 


1 
ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
: 100 
DO.YOU WISH TO SEE THE INTERVENING YEARS? 
1 YES 


0 
100 1 4O7 95 
2 239 2 
3 YT) 10 
4 24 3 
TOTAL 742 307) 113 


DO a Hwee SEE ANY OT RE Rav EAE So 


o OTT PRCGRAM 
1 
Exes ae NUMEER OF THE YEAR YOU WISH TO SEE 
9 00 
Mem OUeWISH TC SEE THE INTERVENING YEARS? 
1 YES” 


0 
900 


WN = 


TOTAL + 


) 692 
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Oe aaa TC SEE ANY OTHER@Y Sano 


1 YES 
v | QUIT PROGRAM 


; 
ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
950 
DO YOU WISH TC SEE THE INTERVENING YEARS? 
: YES | 
0 


a 30 2 
i 


WN = 


53 
32 
11 

uy 
TOTAL 4 2063 


SP ee ee ee ee ee ee ee ee eee 


) 728 


Sh ee 10 SEE ANY O2RER Pharos 


1 VES 
a QUIT PRCGHRAM 


0 
DO. YOU WISH TO SEE THE STEADY STATE VECTOR? 
1 YES 
i 1 


es 
@ 


Pie 
2 


728 


DO YOU WISH TC ENTER DATA? 
0 NC 
7 YES 
° 1 
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BXG ER THE NUMEER OF THE MODEL TYPE 
1 MARKOV HIERARCHICAL 
Z MARKOV LENGTH OF SERVICE 
& MARKOV GENERAL 
4 VACANCY 
5 


ENTER N(INITTIAL STOCK VECTOR) 
300 200 100 
2 Lee eee MATRIX) BY ROWS 
1TH ROW 

oO «15 205 
mvekR 21H ROW 

= (sy, 
ENTER 31H ROW 


uy 
ENTER RPROP (RECRUITMENT PROPORTION VECTOR) 
; wee OF 223 
ENTER ADDITIVE INCREASE 
‘ 25 
ENTER THE PERCENT CODE 
0 NO GRADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 
2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM 


ACh Naa MOU SCITRE TO SEE THE EN TEREDUEATA 2 
NO 
1 Yeo 


I PROPORTION VECTOR 


(ie 


WOULD YQU LIKE 10 CHANGE ANY OF THE DATA? 
5 OTT PROGRAM 
0 
ENTER THE NUMBER OF THE YEAR YOU WISH 10 SEE 
' 10 
DO YOU WISH TC SEE THE INTERVENING YEARS? 
1 YES 


0 
wee pee eeeee ee eee ee eee 
0 1 300 50 
2 200 | " 
3 100 17 
TOTAL 600 100 
10e ed 417 49 
2 88 10 
3 345 47 
TOTAL 850 142 128 


NORA CCe Ra TC SEE ANY OTHER™ YEAR S2 


1 YES 
7 QUIT PRCGRAM 


1 
ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
: 50 


eo Gon ae ale TC SEE THE INTER VENENG YEA RS? 
1 <Y ee 


0 
50 1 827 45 
2 173 3 
3 849 46 
LOTAL 1850 308 Th) Ss 
DO ca a TC SEE ANY OTHER 3 GAnS. 


i hike. PRCGRAM 
1 
ENTFR THE NUMBER OF THE YEAR YOU WISH TO SEE 
100 ; 
DO YOU WISH TC SEE THE INTERVENING YEARS? 
1 YES 
0 
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| 382 


OS Np TOC SEE ANY OTHER YEARS? 


1 WES 
u QUIT PROGHAM 


: 

ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
SU 

WeeYOU WISH TC SEE THE INTERVENING YEARS? 
ye. 


0 
350 


1558 1098 


Cede an mG SEEMANY OTHER YEARS? 
N 


1 YES 
i QUIT PROGRAM 


; 
FNTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
: 400 


ae HWisieO eset DHE INTERVENING YEARS? 


1241 


a aa BO SSEEOANY OTHER RYE Ae Ss 
N 


1 Laas: 
7 ORT PRCG RAH 


0 2 
Oa West 10 SEB THE STEADY STAZE VECTOR? 
NC 
1 Gatos: 
1 


V7 


PERCENTAGES ARE IN STEADY STATE 
999 1 4594 43 
967 9 

3 5039 48 

TOTAL 10600 100) 


EP EP ae SP GE ee Gap “EE ce ce come SRE eS 


ARE YOU THROUGH? 
0 NO 
1 YES 


1 


tes 


1241 


APPENDIX E 
COMPUTER PRINTOUTS OF EXAMPLES FOR OPTION 5 


1. Example (a) 


SIART 
DO YOU WISH TO ENTER DATA? 
0 NC 
1 YES 
: 1 
ENTER THE NUMBER OF THE 4ODEL TYPE 
1 MARKOV HIERARCHICA 
2 MARKOV LENGTH OF SERVICE 
3 MARKOV GENERAL 
4 V ACANCY 
: 1 
ENTER N(INITIAL STOCK VECTOR) 
: 129 74 28 11 
ENTER THE PROMOTION RATE VECTOR. 
THIS VECTOR SHOULD INCLUDE THE PROMOTION RATES FOR THE FIRST 
3 CLASSES. 
: ~102 .046 .033 
ENTER THE WASTAGE PATE VECTOR. 
THe STH VALUE SHOULD INCLUDE THE TOTAL RATE FOR THAT 
STATE DUE TO FITHER PROMOTION OR WASTAGE. 
: Mig 120 <1 .098 
ENTER THE NUMEFR OF THE RECRUIT TYPE 
1 FIXED RECRUIT VECTOR 
2 ALDITIVE (RECRUIT SIZ ) 
3 MULTIPLICATIVE (RECRUT SIZE) 
4 ADDITIVE (SYSTEM SIZE 
5 MULTIPLICATLVE (SYSTEM SIZE) 
: 5 
ENTER RPROP (RECRUITMENT PROPORTION VECTCR) 
: imo00 0 
ENTER MULTIPLICATIVE FACTOR 
; 04 
ENTER THE PERCENT CODE 
0 NO GRADE PERCENTAGES 
1 GRADE SIZE AS PERCENT OF TOTAL GRADE SIZE 
2 GRADE SIZE AS PERCENT OF ORIGINAL GRADE SIZE 
7 QUIT PROGRAM ; 
: 1 
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Eee ace LIKE TO SEE TEE ENTERED SUA 
1 YES 


: 
P MATRIX 
0.728 0.102 0 0 
0 0.83 0.046 0 
0 0 0.867 0.033 
0 0 0 0.902 
N VECTOR 
129 74 28 
OPTION = 
FACT = 
RECRUITMENT PROPORTION VECTOR 
et Oy 0 
WOULD YOU LIKE TC CHANGE ANY OF THE DATA? 
1 YES 
7 QUIT PROGRAM 
: 0 
ENTER THE NUMPER OF THE YEAR YOU WISH IO SEE 
° 10 


ma at a TC SEE THE INTERVENING YEARS 
N 
Ys 


0 
T N PERCENT R 
0 1 129 53 
2 74 31 
3 28 12 
r 17 5 
TOTAL 242 100) 
10 1 148 33) 
2: 82 3 7 
3 Da 10 
4 10 4 
TOTAL 268 111 42 


DO YOU ees TC SEE ANY “OUHERD Yale 
N 


YerS 
a QUIT PKCGRAM 
1 
ENTER THE NUMEER OF THE YEAR  YOUS tS enor 


50 
po vou RISi TO SEE THE INTERVENING YEARS? 
1 YES : 
‘ 0 


830 


50 


Sito 


5 
1 
0 
3) 
2) OZ 
aoe NG” Me@poeeewaANY OTHER YEARS? 


1 Nes) 
i QUIT PRCGEAM 


1 
ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
: 100 
DO YOU WISH TC SEE THE INTERVENING YEARS? 
1 YES 


TOTAL 400 1 


> ee ee cee eee ee ee eee ee eee ee ee eee ee ee ee eee ee 


0 
100 1 362 331 
Z 206 on 
3 66 10) 
= 20 3} 
TOTAL 654 270 101 


MeerOU WISH TO SEE ANY OTHER YEARS? 


0 NC 
1 viES 
7 GUIT PEKCGRAM 
v 
DO YOU VISH TO SEE THE STEADY SIATE VECTOR? 
NC 
1 YES 
' 1 
PERCENTAGES ARE IN STEADY STAIE | 
999 1 362 ( 55) 
2 206 ( “i 
3 06 10 
4 20 3) 
TOTAL 654 100 101 
ARE YOU TYROUGH? 
0 NO 
1 YES 
: 1 
2. Example (b) 
START : 
DO, YOU WISH TC ENTER DATA? 
N 
1 YES 


1 


8 1 


ON MATRIX) BY ROWS 


ENTER N (INITIAL STOCK VECTOR) 


pO oas, 
) 


SEU OI 
i 


J OgE HIE 

Gs 04, a, Se 
BCs Sree 
OV EQ bey A -4 bd FH] 
MD FH UWhI oO 
OZARHEt-! 
Pry bth ibed i! 
Nh Oa -4 
OTmWODtID 


eb-t[2, 6 Fat 


I 
) 
(RECRUITMENT PROPORTION VECTOR) 


S 
REROP 


ag 

[1 

Et em NY i) 
Prat} ee 
fx 


ENTS 5. 


THE. DA TAS 


OF TOTAL GRA 
SZ 


TIVE FACTOR 
OF ORIGINAL 


TO SEE THe ENTERED ee a 


DRO POR?) [ONS VES Lar 


NT 


j 
WOULD YOU LIKE 

NC 

QUIT PRCGKAM 


WOULD YCU LIKE TO CHANGE ANY OF 
YES 


ENTER THE NUMBER OF THE YEAR YOU WISH TO SEE 
10 

Mom. Ou WiSH TC SEE THE INTERVENING YEARS? 
1 Ygs 


0 
oeeeeti(i‘CiéiCOséOCt;CtCN CCR ERCENT OUR 
0 1 300 50) 
2 300 a5) 
3 100 17 
TOTAL 600 100 
10 1 316 48 
2 6a 10 
3 277 13 
TOTAL 24 110 87 


DO YOU WISH TO SEE ANY OTHER YEARS? 
0 NC 
1 YES 

7 QUIT PRCGRAM - 

: 1 

ENTER THE NUMEER OF THE YEAR YOU WISH TO SEE 
i 50 

DO YOU WISH TO SEE THE INTERVENING YEARS? 

N 

1 YES 


0 
59 


lal) ana 


TOTAL 


P23 


OS eae TE SEE ANY OTHER YEARS? 


1 CaS) 
Ty QUIT PROGRAM 


FNTER THE NUMEES CF THE YEAR YOU VISH TO SEE 
: 100 ‘ 
ar no TO SEE THE INTERVENING YEARS? 
N 


ey ES 


0 
100 


2O.3 


&3 


DO YOU one TO SEE ANY OTHER YEARS? 


= 


ES 
gULT PRCGRAM 


~~) 2 O 


DO YOU WISH TC SEE THE STEAD SrA Tee ele. 
1 YES 
1 
PERCENTAGES ARE IN S 
=o) oh) 1 


Z0S 


a SS Se SS Se ce ee ee ee ee eee ee eee Se eee ee 


ARE YOU THROUGH? 
0) NC 
1 YES 
1 
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APPENDIX F 


APL COMPUTER FROGRAM FOR STEADY STATE DISTRIBUTION FUNCTION 


This 1s the APL function that is used to find the steady 


State dastributicn vector for all options. 


Y STDST;Q; IDENT; RV;W; NSS; NHAT1;NHAT2; FACTI ; IQ; ZERO, 
<N 4 


ave (4, ERE BCG, 

Derek) 0 roe 

W<(1,K)o1- (+ P 
) 
P 


>(TYPE= u 5 
FIX:N<R+. pe 
NSS<No Gey N ) 
N~0 ROUND N 


SEND 
ADDREC : NHAT1<RPROP+. x8 (IDENT- P) 
NSS<NHATI3 (+/NHAT1) 


UPR PACLI<CFACTXI DENT 
NHAT2<R+.xf(FACTI-P 
Roe Naal 2% (+/NHAT2 ) 


aa «kot Q 


aleaade inci i aaa S2S 


NSS«ZER +.xBIQ 
END 


MULTSYS FACTI<FACTxIDENT 
NSS<(FACT-1)xRPROP+.x3( FACTI-Q) 
END: TCOUNT<=399 
PERCENT <I (0. 5+100xNSS) 
POTPERCENT<100 
OUTPUT 
— 


MNMNMNNMDNNNNNRP PRP RPPRPRPRPRPRPOONMUF WHR 


WOONONEWNHOROWONMNUE WHF O— 
hy 
$O 
a 
re 
0D 
me 
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